Robust strategies for lossy quantum interferometry 
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We give a simple multiround strategy that permits to beat the shot noise limit when performing 
interferometric measurements even in the presence of loss. In terms of the average photon number 
employed, our procedure can achieve twice the sensitivity of conventional interferometric ones in 
the noiseless case. In addition, it is more precise than the (recently proposed) optimal two-mode 
strategy even in the presence of loss. 
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PACS numbers: 03.65.Ta,06.20.Dk,42.50.St 

The shot noise limit is the minimum noise level that 
the Heisenberg uncertainty relations permit to achieve 
when classical states are employed in the apparatuses. 
Many quantum strategies have been proposed to beat 
the shot noise [H, [3| and to achieve the ultimate Heisen- 
berg limit , but virtually all of them are very sensitive 
to noise and loss of photons f3| . Only very recently some 
interferometric strategies were presented that can beat 
the shot noise even in the presence of relevant losses of 
photons [5, 01 • These are all instances of parallel strate- 
gics [2] , where both arms of the interferometer are sam- 
pled at the same time using a mode-entangled quantum 
state of the light (see Fig. [T^). In addition to the par- 
allel strategies, in quantum metrology it is also possible 
to achieve the Heisenberg limit using multiround (or se- 
quential) strategies 0, 0] which, in the noiseless case, are 
equivalent in terms of resources and of achievable preci- 
sion [§]. 



o' 

C3 



9 



3 



3 

3 



prepa- 
ration 



measu- 
rement 



9 



U 



a) 



b) 



FIG. 1: a) Parallel strategies for interferometry. Both arms 
of the interferometer are sampled at the same time using a 
two-mode entangled state. The phase factor ip is imprinted in 
the state as a phase difference between the two modes, b) Se- 
quential strategy proposed here. A loss-resistant single mode 
state is sent through the first interferometer arm, sampling 
the phase ip. Then a unitary transformation is applied and 
the state is sent back through the other interferometer arm. 
This is needed so that the final phase shift experienced by 
the state is only the relative phase in the interferometer. The 
state is measured after one (or more) round trips. 

Here we detail how multiround strategies can be used 
to perform interferometry — see Fig. [TJd. An appropriate 
input state is prepared (we will analyze two examples be- 
low). This state is fed into the first interferometer arm. 
It picks up a phase f + i!), where is the interferometric 
phase we want to estimate, and ■& is the absolute phase 



picked up by the free evolution in the arm (which is equal 
to the phase which would be picked up also in the refer- 
ence arm) . The main trick of multiround interferometry 
is the use of the unitary 
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where \n) is the Fock basis and M is the largest nonzero 
component of the initial input state. The purpose of this 
unitary is to permute the first M components of the Fock 
state expansion of a state in such a way that, when the 
state is sent back through the reference arm, the absolute 
phase {} is removed from the state (only an irrelevant 
global phase factor e**^'' persists). Thus, at the end of 
the round trip of Fig. [TJs (multiple round trips are also 
possible), only the relative phase ip is imprinted on the 
state. A measurement is finally performed to estimate 
this phase. 

This multiround interferometry employs the same av- 
erage energy and the same modes as the conventional 
(parallel) strategies, but it can achieve twice the sensi- 
tivity and it is more robust against noise. In fact, we 
show that, with an appropriate choice of inputs, our pro- 
tocol permits to estimate the phase with an error which 
is smaller than what is achieved by the strategy detailed 
in Ref. [3] , which is claimed to be the optimal two- mode 
strategy. In the presence of loss, this is not unexpected, 
as parallel protocols rely on entanglement which is no- 
toriously fragile to noise. Instead, in the noiseless case, 
this can be seen easily with a simple example. Recall 
that the Heisenberg limit I4I is essentially an applica- 
tion of the time-energy uncertainty, A(/3A/i ^ 1/2, where 
h is the generator of the unitary that inserts the phase 
if into the system IJ], namely h = a) a (a being the 



annihilation operator of the first arm of the interferome- 
ter). Optimal two-mode states, such as the NOON state 
(|A^0) + \m))/^/2, have A/i = N/2. However, the corre- 
sponding single-mode "NO" state of same average num- 
ber of photons, namely the state (|27V) -I- |0))/'\/2, has 
Ah = N . Both states achieve the Heisenberg-limited 
sensitivity of = l/(2Aft,), but in terms of the aver- 
age number of photons N , the NOON state can achieve 
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'^'Pnoon = 1/-^, whereas the NO state can achieve 
^fNO = l/(2-^V), i.e. twice the sensitivity. The NO state 
is, however, just as sensitive to noise as the NOON state: 
the loss of a single photon renders both states useless to 
phase estimation. 

The rest of the paper is devoted to presenting two ex- 
amples of loss-resistant multiround interferometry based 
on two different input states, the optimal phase state and 
the single-mode M&M state. Both are able to beat the 
optimal parallel strategy for some values of parameters. 

We start by considering the single-mode optimal phase 
state l^opt) introduced in [^Q (and later extended to the 
multi-mode case in [loj, building on fll']), i.e. the state 
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where M is a parameter identifying the average num- 
ber of photons N — M/2. In the noiseless case, this 



state can be used to achieve the ultimate precision [12 1 
in the absolute phase estimation, i.e. the Heisenberg scal- 
ing ~ 1/N. Moreover, this state is highly robust to 
loss, since the loss tends to deplete primarily the com- 
ponents with large Fock numbers n which are not very 
populated in this state. Following the scheme of Fig. [T]d, 
this state is sent through the first arm of the interfer- 
ometer, where it is subject to a phase shift of (f + 'd, 
where ip is the interferometric phase we want to esti- 
mate, and "i? is the absolute phase picked up by the free 
evolution in the arm. During this transit, the state is also 
evolved by the loss map, described by the Kraus opera- 
tors Ki = (77^^ — l)*/^a''77°*°/^/\/iT, where a is the annihi- 
lation operator of the optical mode and 77 is its transmis- 
sivity or quantum efficiency. (Note that the loss and the 
phase accumulation commute, so that the order in which 
we apply these two transformations is irrelevant.) Then, 
the state is subject to the unitary evolution U of Eq. ((J) 
and it is sent back along the reference arm. Thanks to the 
permutation of the Fock components that U applies to 
the state, the free evolution is effectively reversed, so that 
the absolute phase that was picked up in the first arm 
is removed while the radiation travels back through the 
reference arm. Also in the reference arm the state is typ- 
ically subject to the loss (although there are interesting 
cases where the reference may be considered noiseless). 
Finally, the state is subject to the measurement. For the 
optimal phase state, a good measurement Q can be ob- 
tained by considering the orthogonal POVM composed 
by the projectors on the Pegg-Barnett states 
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with I = 0, • • • , M. The RMS of the probability dis- 
tribution obtained from this POVM is a function of the 
interferometer phase (p. Its minimum value gives the min- 
imum error that our interferometer can achieve, which is 



plotted as the continuous lines in Fig. [2l In addition, we 
have also directly estimated the error through the Holevo 
variance [3 A$ = (5^^ - 1)3, where 5$ = |(e'*)| is the 
average value of the function e** weighted with the prob- 
ability obtained from the continuous POVM |$)($| d'^, 
obtained using the phase states of Eq. ^ for arbitrary $. 
The Holevo variance is more appropriate than the RMS 
for the estimation of the phase error, since the phase is 
a periodic quantity [l3| . However, as is clear from our 
plots, the Holevo variance is well approximated by the 
RMS when these two quantities are small enough (com- 
pared to 2it). 

It is tedious but straightforward to calculate that, af- 
ter a round trip characterized by the same quantum effi- 
ciency 77 in both arms, the state IV'opt) is evolved into 
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where the second sum runs between max(0, i — j) and 
M — j. From this state, one can estimate the probability 
distribution of the POVM of Eq. namely 
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The error in the estimation of the phase from a measure- 
ment on p is given by the RMS oip{l), plotted in Fig.[2]as 
a function of the average photon number N and quantum 
efficiency 77. 

The second input state we consider is a single-mode 
analogous of the M&M state introduced in [H], i.e. the 

state 
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where M > M' and whose average photon number is 
N = {M + M')/2. Again, it is straightforward to obtain 
the output state, after the round trip of Fig. [TJa: 

^ - Ef=-s^j\j + s){j + s\ + E'Lom{j\ (5) 
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with 5 = M - M' imd fij = (1 - ,y)2i-j^A/-i+j. To 
extract the phase from the state a, in analogy to the 
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FIG. 2: Left) Solid line: error in the estimation of the phase 
from the state p, i.e. minimum (over ip) of the RMS of the 
probability p(/), as a function of the average number of pho- 
tons A'^ for 77 = 0.9. Lower dashed line: plot of the error 
from the optimal two-mode lossy interferometry from Q . Our 
single-mode method achieves a higher precision for a wide 
range of parameters. Plus signs: Holevo variance A$, which 
closely tracks the minimum RMS. Dot-dash line: behavior of 
the two mode state (|A''0) + |07V)/%/2, which is optimal only 
for high r] and low A'^. Upper grey dashed line: average over 
tp of the RMS. The gray shaded area encloses the "quantum" 
trajectories, i.e. the ones included between the shot noise limit 
1/y/Nri and the Heisenberg limit 1/A'^. Right) Same as the 
previous, but as a function of ri for TV = 20. 



two- mode case of |5[, we can measure the observable 

M' 
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fc=0 

Then, the error on the phase can be obtained from the 
RMS of A using error propagation, namely 
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where 9 = J2kLo ak+ak-5+Pk+f3k+5 and T = X^feio 
This quantity is plotted in Fig. [3] from which it is evi- 
dent that, also in this case, the multiround protocol can 
achieve a better sensitivity than the optimal two-mode 
one. The single-mode M&M state can take advantage of 
most of the implementation ideas presented in Q for the 
two-mode case. It may appear surprising that the state 
a permits to achieve a greater precision than the NOON 
state i\NO) + \0N))/V2 even for high values of 77. As 
discussed above, this is essentially due to the fact that a 
single mode state performs better than a two-mode state 
in terms of the average number of photons N. 

In conclusion, we have given a strategy for determining 
the relative phase in an interferometer using single-mode 
states that are sent through the interferometer in a round 
trip, interleaved by the unitary U of Eq. ([1]). This entails 
that a) in the noiseless case a double sensitivity can be 
reached over the optimal two-mode states (such as the 
NOON state) in terms of the average number of pho- 
tons, b) in the lossy case we can achieve a better phase 
sensitivity than what is claimed to be the optimal two- 
mode strategy fij , proving that multiround protocols are 
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FIG. 3: Left) Solid line: error in the estimation of the phase 
from the state a of Eq. ([5}, i.e. minimum (over ip) of the 
function Aip of Eq. pop . as a function of the average number 
of photons iV for 7? = 0.9. Here M' = 3 and Af = 2iV - 3. 
Dashed line: plot of the error from the optimal two-mode lossy 
interferometry from [J . Again, as in Fig. [2J our single-mode 
method achieves a better sensitivity in some range. Dot-dash 
line: behavior of the state (|iV0) + \ON)/y^. The inset is 
an enlargement for small values of N. Right) Same as the 
previous, but as a function of -q for N = 20, M = 30, M' — 10. 



preferable in the presence of noise. The robustness in 
the face of loss stems from two main properties. On one 
side there is no entanglement between different modes, 
and it is well known that entanglement is very sensitive 
to noise [3]. On the other side, the fact that we are us- 
ing a single mode permits to double the phase sensitivity 
over the two-mode entangled case, since all the photons 
travel through one mode only. One may object that the 
increased phase sensitivity arises because we are devoting 
more resources to the estimation. This objection is un- 
founded since the average number of photons employed 
in the two strategies is the same, N. One cannot even 
say that in the two-mode strategies the phase ip is sam- 
pled by less photons, because the number of photons that 
travel through an arm of an interferometer is an unde- 
fined quantity (the "which path" information is comple- 
mentary to the phase information). One can only bound 
the number of photons traveling through each arm with 
the total number of photons, N , injected in the interfer- 
ometer. 

We thank R. Demkowicz-Dobrzanski for having kindly 
provided the data of the optimal two- mode state of [3| ■ 
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